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Extensively applied to both light and heavy meson decay and standing as one of the most successful 
strong decay models is the ^Po model, in which qq pair production is the dominant mechanism. 
The pair production can be obtained from the non-relativistic limit of a microscopic interaction 
Hamiltonian involving Dirac quark fields. The evaluation of the decay amplitude can be performed 
by a diagrammatic technique for drawing quark lines. In this paper we use an alternative approach 
which consists in a mapping technique, the Fock-Tani formalism, in order to obtain an effective 
Hamiltonian starting from same microscopic interaction. An additional effect is manifest in this 
formalism associated to the extended nature of mesons: bound-state corrections. A corrected ^Po 
is obtained and applied, as an example, to 6i ujn and ai pn decays. 
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I. INTRODUCTION 

A great variety of quark-based models are known that 
describe with reasonable success single-hadron proper- 
ties. A natural question that arises is to what extent a 
model which gives a good description of hadron prop- 
erties is, at the same time, able to describe the com- 
plex hadron-hadron interaction or by the same principles 
hadron decay. In particular, the theoretical aspects of 
strong decay have been challenged by QCD exotica (glue- 
balls and hybrids) where a consistent understanding of 
the mixing schemes for these states is still an open ques- 
tion 3 1 . The nature of the family of "new mesons" 
X,Y,Z 4| is another unsolved puzzle: are they actually 
new qq mesons, hadronic molecules or something else? In 
the direction of clarifying these questions is the successful 
decay model, the "^Pq model, which considers only OZI- 
allowed strong-interaction decays. This model was intro- 
duced over thirty years ago by Micu [H and applied to 
meson decays in the 1970 by LeYaouanc et al. Q. This 
description is a natural consequence of the constituent 
quark model scenario of hadronic states. 

T. Barnes et al. have made an extensive sur- 

vey of meson states in the light of the ^Pq model. Two 
basic parameters of their formulation are 7 (the inter- 
action strength) and f3 (the wave function's extension 
parameter). Although they found the optimum values 
near 7 = 0.5 and (3 = 0.4 GeV, for hght IS and IP de- 
cays, these values lead to overestimates of the widths of 
higher-L states. In this perspective a modified qq pair- 
creation interaction, with 7 — 0.4 was preferred. 

In the present work, we employ a mapping technique 
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in order to obtain an effective interaction for meson de- 
cay. A particular mapping technique long used in atomic 
physics [llj, the Fock-Tani formalism (FTf), has been 
adapted, in previous publications p^-[T6^, in order to 
describe hadron-hadron scattering interactions with con- 
stituent interchange. Now this technique has been ex- 
tended in order to include meson decay. We start from 
the microscopic qq pair-creation interaction, as will be 
shown, in lower order, the ^Pq results are reproduced. 
An additional and interesting feature appears in higher 
orders of the formalism: corrections due to the bound- 
state nature of the mesons and a natural modification in 
the qq interaction strength. 

In the Fock-Tani formalism one starts with the Fock 
representation of the system using field operators of 
elementary constituents which satisfy canonical (anti) 
commutation relations. Composite-particle field opera- 
tors are linear combinations of the elementary-particle 
operators and do not generally satisfy canonical (anti) 
commutation relations. "Ideal" field operators acting on 
an enlarged Fock space are then introduced in close cor- 
respondence with the composite ones. Next, a given uni- 
tary transformation, which transforms the single compos- 
ite states into single ideal states, is introduced. Applica- 
tion of the unitary operator on the microscopic Hamilto- 
nian, or on other hermitian operators expressed in terms 
of the elementary constituent field operators, gives equiv- 
alent operators which contain the ideal field operators. 
The effective Hamiltonian in the new representation has 
a clear physical interpretation in terms of the processes 
it describes. Since all field operators in the new repre- 
sentation satisfy canonical (anti)commutation relations, 
the standard methods of quantum field theory can then 
be readily applied. 

In this paper we shall extend the FTf to meson de- 
cay processes. In the next section we review the basic 
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aspects of the formalism. Section Hill is dedicated to ob- 
tain an effective decay Hamiltonian. In section HVl two 
light mesons decays examples are calculated bi — > um 
and fli — > pn. The summary and conclusions are fol- 
lowed by appendixes which detail the method employed 
throughout this work. 



II. MAPPING OF MESONS 

This section reviews the formal aspects of the mapping 
procedure and how it is implemented to quark-antiquark 
meson states [l^. The starting point of the Fock-Tani 
formalism is the definition of single composite bound 
states. We write a single-meson state in terms of a meson 
creation operator as 



Ml\0), 



(1) 



where |0) is the vacuum state. The meson creation op- 
erator is written in terms of constituent quark and 
antiquark creation operators and , 



(2) 



^i^" is the meson wave function and g^|0) — qi,\0) — 
0. The index a identifies the meson quantum numbers 
of space, spin and isospin. The indices and i' denote 
the spatial, spin, flavor, and color quantum numbers of 
the constituent quarks. A sum over repeated indices is 
implied. It is convenient to work with orthonormalized 
amplitudes. 



(3) 



The quark and antiquark operators satisfy canonical an- 
ticommutation relations, 

i-^} = = Ut^, = {qti,ql} = o. (4) 

Using these quark anticommutation relations, and the 
normalization condition of Eq. ([3]), it is easily shown that 
the meson operators satisfy the following non-canonical 
commutation relations 



where 



In addition, 



[g^,Mt] = ci>^-9-t , 



[g-.,Mt] = -<i>^J-gt 



(6) 



(7) 



The presence of the operator I^afj in Eq. ([5]) is due 
to the composite nature of the mesons. This term enor- 
mously complicates the mathematical description of pro- 
cesses that involve the hadron and quark degrees of free- 
dom. The usual field theoretic techniques used in many- 
body problems, such as the Green's functions method. 



Wick's theorem, etc, apply to creation and annihilation 
operators that satisfy canonical relations. Similarly, the 
non-vanishing of the commutators [g^,Mj^] and [gi,,M^] 
is a manifestation of the lack of kinematic independence 
of the meson operator from the quark and antiquark op- 
erators. Therefore, the meson operators Ma and Mj^ are 
not convenient dynamical variables to be used. 

A transformation is defined such that a single-meson 
state I a) is redescribed by an ("ideal") elementary- meson 
state by 



<|0), 



(8) 



where an ideal meson creation operator. The ideal 
meson operators toJ, and TOq satisfy, by definition, canon- 
ical commutation relations 



(9) 



The state |0) is the vacuum of both q and m degrees of 
freedom in the new representation. In addition, in the 
new representation the quark and antiquark operators 
q^ , q, q^ and q are kinematically independent of the mj, 
and TUa 

[q^,ma] = [g^, mjj] = [g^, nia] = [g^, ml] = 0. (10) 

The unitary operator U of the transformation is 

C/(i) =exp[tF], (11) 

where F is the generator of the transformation and t a 
parameter which is set to —tt/2 to implement the map- 
ping. The next step is to obtain the transformed oper- 
ators in the new representation. The basic operators of 
the model are expressed in terms of the quark operators. 
Therefore, in order to obtain the operators in the new 
representation, one writes 



q{t)^U-^qU, q{t) = U-^qU. 

The generator F of the transformation is 



[Ma,M^] = 0, (5) where 



with 



F^ml Ma - Mir 



1=0 

[Ma,Ml]^Sa0 + 0($ 

[Ma,Mf3]^[MlMl]=0 



(12) 



(13) 



(14) 



n+l^ 



(15) 



It is easy to see from (|13l) that F^ ~ —F which ensures 
that U is unitary. The index i in (|14p represents the 
order of the expansion in powers of the wave function 
<i>. The Ma operator is determined up to a specific order 
n consistent with (jlSp . The examples studied in [l2| 
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required the determination of Mq*^ up to order 3 as shown 
below 



Mi") 



Mi' 







M(2) = ^A^pMf, ; M(3) ^ T„;5^ M^, (16) 

withTc^^ = — [Ma, A^^]. In the "zero-order" approx- 
imation, overlap among mesons is neglected and terms 
of the same power in the bound-state wave function 
i^a) collected. In order to have a consistent 
power counting scheme, the implicit entering 
via Eq. ([2]) are not counted. The consequence of this is 
that the equations for TOq and are manifestly sym- 
metric, 



drUa (t) 

Jt 
dMait) 

Jt 



= [Mait),F]=-ma{t), (17) 



and their solutions involve only trigonometric functions 
oft, 

ITT-ait) = MaSint + maCOSt, 

Ma{t) = MaCOSt — rriaSuvt. (18) 

The equations of motion for the quark operators q and 
q can be obtained by making use of Eq. ([7]) in a similar 
way. 



dqf,{t) 
dt 



dqt,{t) 
dt 



= [q^{t),F].m 



In the zero-order approximation, the effects of the meson 
structure are neglected resulting 



g(«)(i) = gp, e\t)=q., 

'^a ^ (*) — ™a COS t + Ma siu t, 

M^°) {t) = Ma cos t-TTia sin t. 
In first order one has 



(20) 



= - ["la sint + Ma{l~ cost)] , 

qij-){t) = ^^I'^qt^ [ma sint Af„ (1 - cost)] , 
™W(t)=0, MW(t) = 0. (21) 

The second and third order solutions to were calcu- 
lated in reference T5| and appear again, for completeness, 
in appendix [XI together with the higher order operators 
required in our calculation. 

Once a microscopic interaction Hamiltonian H is de- 
fined, at the quark level, a new transformed Hamiltonian 
can be obtained. This effective interaction we shall call 
the Fock-Tani Hamiltonian and is evaluated by the ap- 
plication of the unitary operator U on the microscopic 
Hamiltonian TCft = U~'HU . The transformed Hamil- 
tonian TipT describes all possible processes involving 



mesons and quarks. The general structure of T^ft is of 
the following form 



TYft ~'H„+ Tim + Tir 



(22) 



where the first term involves only quark operators, the 
second one involves only ideal meson operators, and Timq 
involves quark and meson operators. 

In TipT there are higher order terms that provide 
bound-state corrections (also called orthogonality cor- 
rections) to the lower order ones. The basic quantity 
for these corrections is the bound-state kernel A{pT; Xv) 
defined as 



A{pT;Xu)^<pPJ<P*^ 



(23) 



To discuss the physical meaning of the bound-state 
corrections and how they modify the fundamental quark 
interaction we shall present an example, in a toy model 
similar to the model studied in [l3], where the basic ar- 
guments are outlined. In this example, the starting point 
is a two-body microscopic quark-antiquark Hamiltonian 
of the form 



H2q = T{fi)qlqf, + T{i^)qlq,, + Vqq{fii^;(Tp)qlqlqpq, 



(24) 

The transformation TipT = U^' H2qU is implemented 
again by transforming each quark and antiquark operator 
in Eq. ([24]) . where a similar structure to Eq. ([22]) is 
obtained. In free space, the wave function $ of Eq. 
satisfy the following equation 

77(/i^.;ap)$-P = e[,]$q, (25) 

where H^fiiy; ap) is the Hamiltonian matrix 

H{^„.■,ap) = S^[,]S,[,][Ti[a])+T{[p])] 



+Vqq{pi^;ap), 



(26) 



£[0,] is the total energy of the meson. There is no sum 
over repeated indices inside square brackets. 

The effective quark Hamiltonian Tiq has an iden- 
tical structure to the microscopic quark Hamiltonian, 
Eq. ([M]) . except that the term corresponding to the 
quark-antiquark interaction is modified as follows 



Va 



[Vqq-H A-AH + AH A] 



(27) 



where Vqq = Vqq{pv]ap) and the contraction H A = 
H{ph';T^) A{T^;ap). An important property of the 
bound-state kernel is 



Aip,y;ap)<P^a'^K'', 



(28) 



which follows from the wave function's orthonormaliza- 
tion, Eq. Q. In the case that <I> is a solution of Eq. (|25p . 
the new quark-antiquark interaction term becomes 

V,,{p,.; ap) = V,,{p,y; ap) - ^ ea<i>TK'- (29) 
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The spectrum of the modified quark Hamiltonian, Tin, is 
positive semi-definite and hence has no bound-states . 
This result is exactly the same as in Weinberg's quasipar- 
ticle method , where the bound-states are redescribed 
by ideal particles. The new Vqq is a weaker potential, 
modified in such a way that no new bound-states are 
formed. 

In the quark- meson sector of Eq. (|22p in Timq appears 
a term related to spontaneous meson break-up 



'm^qq 



(30) 



with 



(31) 

Again, in the case that <i> is a solution of Eq. a 
straightforward calculation demonstrates that H^^qq — 
0. When there is no external interaction, this result is a 
direct consequence of the bound-state's stability against 
spontaneous break-up. This term can be interesting in 
studies related to dense hadronic mediums. For these 
systems the wave function is, in general, not a solution 
of Eq. (pSj) and the strength of the potential V{p,v; a) is 
now only decreased [Tsj . 

In the ideal meson sector Tim many similar approaches 
to FTf [l^ have obtained the meson-meson scattering in- 
teraction in the Born approximation: Resonating Group 
Method (RGM) [l^. Quark Born Diagram Formalism 
(QBDF) 



(32) 



where Tmm is the kinetic term and Vmm is the meson- 
meson interaction potential with constituent interchange. 
This potential is given by 



dir 



(33) 



where V^*^ is the direct potential (no quark interchange). 



T/'exc 
^ mm 



the quark exchange term and the intra- 



exchange term. As shown in Ref. [fj] and [13|, if one 
extends the FT calculation to higher orders a new meson- 
meson Hamiltonian is obtained 



5H„ 



(34) 



where SHm„i is the bound-state correction Hamiltonian. 
If the wave function $ is chosen to be an eigenstate of the 
microscopic quark Hamiltonian, then the intra-exchange 
term V"^* is cancelled 



SHm.rn. — 0. 



(35) 



In summary, these examples reveal an important and 
common feature of bound-state corrections: they weaken 
the quark-antiquark potential. In the next section we 
shall follow the same procedure for a quark pair creation 
interaction, which is fundamental for the description of 
meson decay. Similar to the toy model, the resulting 
interaction that describes meson decay, will contain a 
Born order contribution and a bound-state correction. 



III. 



THE ^Po DECAY MODEL IN THE 
FOCK-TANI FORMALISM 



In the paper of E. S. Ackleh, T. Barnes and E. S. Swan- 
son 01 a formulation of the '^Pq model is presented. It 
regards the decay of an initial state meson in the pres- 
ence of a qq pair created from the vacuum. The pair pro- 
duction is obtained from the non-rclativistic limit of the 
interaction Hamiltonian Hj involving Dirac quark fields 



(36) 



where 7 is the pair production strength. For a qq meson 
A to decay to mesons i? + C we must have {qq)A 
{qq)B + {qq)c- To determine the decay rate a matrix 
element of ([55]) is evaluated 



{BC\Hi\A) = 5{Pa -Pb- Pc) hf^. 



(37) 



The evaluation of h fi is performed by diagrammatic tech- 
nique for drawing quark lines. The hfi decay amplitude 
is combined with relativistic phase space, resulting in the 
differential decay rate 



dV 



A^BC n_r,^BEc,, |2 



= 2t:P- 



Ma 



-\hf. 



(38) 



which after integration in the solid angle f2 a usual choice 
for the meson momenta is made: Pa = (P = \Pb\ = 
\Pc\)- 

In our approach, the starting point for the Fock-Tani 
hfi is also the microscopic Hamiltonian Hi in (|36p . The 
momentum expansion of the quark fields, color and flavor 
are not represented explicitly, is 



(27r)3 

+v,{-k)qj{-k)]e''^-^. 



(39) 



In the product ip{x)ip{x) we shall retain only the q'^q'' 
term, which yields from Eq. (j36[) a Hamiltonian in a 
compact form. 



Hi = V^,^ qlql 



(40) 



where sum (integration) is again implied over repeated 
indexes. In the compact notation, the quark and anti- 
quark momentum, spin, flavor and color are written as 
A* = iPfi,Sfj,, ff_t,Cf_,); V = ip„,s^, f^,c^), while the pair 
creation potential Vf^i, is given by 

Vf,^ = 2m,7 S{pf,+p,,)u,^f^c^{pf,)vs^f^cAPi^)- (41) 



It should be noted that since Eq. p6|) is meant to be 
taken in the nonrelativistic limit, Eq. (I41|l should be as 
well. In the meson decay calculations, of the next section, 
this limit is considered. 
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Applying the Fock-Tani transformation to Hi one ob- 
tains the effective Hamiltonian 



1-Lft = Hi U. 



(42) 



The physical quantities in the FTf appear in a second 
quantization notation. The effective decay amplitude will 
be a product of the ideal meson operators with the fol- 
lowing structure in the ideal meson sector: m^ra^ra. To 
obtain this product corresponds to expand, in powers 
of the wave function, up to third order. A Hamiltonian 
that describes this decay process, which we shall call H^, 
can be extracted from the mapping by the following 
products 



H — V n^^'^'nT 



t(3)^t(0) 



(43) 



After the substitution of Eqs. jH]), jj!]) and 

into ([331) results in the effective meson decay Hamiltonian 



H,n = /^"(a, 7) V^^ ml^m''pmj 



where 



Q /3 7 • 



(44) 



(45) 



In the ideal meson space the new initial and final states 
involve only ideal meson operators |^) = to^^IO) and 

\BC) = mj,mt|0). The ^Pq amplitude is obtained in 



the FTf by an expression equivalent to Eq. (|37|) , 

{BC\Hi\A) = (0|m,m/3i/„mt|0) 

= r'^(a,/3,7)V + /'''(/3,a,7)^M- (46) 

The term /^''(/J, 0,7) of (gH) is shown in Fig. ([TJi), the 
term /^"(a, /3, 7) corresponds to the same diagram with 
a ^ /?. 

In the FTf perspective a new aspect is introduced to 
meson decay: bound-state corrections. The lowest order 
correction is one that involves only one bound-state ker- 
nel A(^v,ap). This implies that the Hamiltonian repre- 
senting this correction must be of fifth order in the power 
expansion of the wave function. 

We shall call this new Hamiltonian, with the same ba- 
sic operatorial structure ml^mpirij, of 6Hm- The only 
combinations that results in a fifth order Hamil- 

tonian are 



SH„i — 



-e^q'J'^ + S'^q'J''^ + q^^r^] V,.. (47) 



Details of this calculation is found in the appendix IB] 
The bound-state corrected ^Pq Hamiltonian, which shall 
be called the C^Pq Hamiltonian, is 



H 



C3P0 



Hm + SHm 



(48) 



where yCSPO jg condensed notation for 

1 



V 



C3P0 



+ A(pt;/i(t) 



err 



(49) 



IV. LIGHT MESON DECAY EXAMPLES 

The light meson sector is an interesting test ground 
where the effects of the bound-state correction can be 
compared to the usual ^Pq model. In particular, as 
examples, two specific decay processes will be studied: 
bi — !■ con and ai pn. The wave function and details of 
the matrix elements are found in the appendix [Cl The 
general decay amplitude can be written as 



lC3P0 



1 >lC3P0 



(50) 



For the first decay process, bi — > cjtt, results in a decay 
amplitude given by 



A1C3P0 



with 



Coi = — 
C21 = - 



35/2 



1 - -x' 
9 



ei{x) 



75/23 



21 



2I1/2 



39/2 

where x = P/ [3 and 



eiix) 



217/2 

7V29 



62(2:) 



61(2::) = exp 



X 

'V2 



62(2;) — exp 



%x^ 
' 28 



(51) 

e2{x) 
(52) 

.(53) 



The decay rate has a straightforward evaluation, by sub- 
stituting (ini), ([52]) in ((501) and then in obtaining 



- h\ — yUJTT 



2V^x^j' (C^oi+Cl,). (54) 



Mb, 



The second decay process, oi — > pir, is similar to the 
former one and results in the following amplitude 



M 



C3P0 



with 



Coi — 



29/2 

3572 

n — 2 
C21 = ~X 



1 - -X^ 

9 

39/ 



Col>oo(^^.) + C2lr2o(^^.;) , (55) 

211/2 



ei(x) 



75/23 



1 

21^ 



2^5 



2ei(^)-^j2^e2(a;) 



32 VI 

and by a similar procedure one obtains 



2V7r: 



EnE 



M„ 



7' (C^oi+Cl) 



e2{x) 
(56) 

(57) 
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FIG. 1: Diagrams representing the C^Po model. Diagram (a) corresponds to the ^Po amplitude. Diagrams (b), (c) and (d) are 
the bound state corrections. The complete hfi amplitude includes the diagrams above plus diagrams with q ^ /3. 



In the former equations, 62 (x) — 0, recovers the original 
•^Pq results. 

In addition to the decay widths F, bi and oi mesons 
have D/S ratios, which give a sensitive test for decay 
models. By definition, these quantities are obtained from 
the ratios of C21 and Coi coefficients, in equations ((52|) 
and dnH). 



D 




9 

— X 


"2I/2 , , 01/235/2 r , ■ 

^ei{x) ~ ^ ^e2(x) 






ai— >p7r 


[1 - 
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,3/2 a^/^s^/^ / \ 

^ei{x) - 62(3;) 



[1 



ei(x) 



(58) 



75/2 



The meson masses assumed in the numerical calculation 
were M„ = 138 MeV; Mp = 775 MeV; M^, = 1230 MeV; 
Mb, = 1229 MeV; = 782 MeV 0. 

The choice of SHO wave functions allow exact eval- 
uations of the decay amplitudes even in the corrected 
model. A first new aspect that appears is the presence 
of a new dependence in the exponential of the corrected 
term. This implies in a different range for the bound- 
state correction due to the fact that e2{x) / ei{x) ^ as 
X — > 00. 

The correction introduces the bound-state kernel, Eq. 
to the calculation of the decay processes. A gen- 



eral sum over the meson index a is present and as stated 
before, this index represents the quantum numbers of 
space, spin and isospin. The OZI-allowed decays rep- 
resent, flavor conserved continuous (anti)quark lines. A 
direct consequence of this fact is the possibility to sum 
over a larger set of mesons in the a index. In our cal- 
culation the sum was restricted only to the final state 
mesons. In the — > ujtt'^ decay, there are two bound- 
state kernel contributions one associated to uj meson and 
the other to 7r+. Similarly, the — > /9+7r° decay has 
two bound-state kernel contributions one associated to 
/9+ meson and the other to the 7r°. 

In this example, the parameters were chosen in order 
to give a closer fit to the experimental data. In the 61 
decay, width and partial waves are known with accuracy. 
The "^Po model's optimum fit for the bi data (F and D/S 
ratio) is achieved with 7 = 0.506 and (3 = 0.397 GeV. In 
the C^Pq rnodel a similar fit is obtained with 7 ~ 0.539 
and P = 0.396 GeV. These parameters are used in the 
two models to describe the ai decay. The results for F 
as a function of (3 appear in figure [2] and specific values 
are presented in table [D In figure [3l the D/S ratios for 
the two models are plotted. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have presented an alternative ap- 
proach for meson decay which consists in a mapping 
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TABLE I: Decay rates ^Po (7 = 0.506 e /J = 0.397 GeV ) and C^Po (7 = 0.539 e /3 = 0.396 GeV ) 







r (MeV) 






D/S 




Decay 


Exp [20] 


'Po 




Exp [20] 


'Po 


C^Po 


bi —* UJTV 


142 


143 


142 


0.277(27) 


0.288 


0.273 


ai piv 


250 to 600 


543 


543 


-0.108(16) 


-0.149 


-0.113 





FIG. 2: Decay rates for 61 ujtt and ai pn decays, for Po 
(7 = 0.506) and C^Po (7 = 0.539) models . 



FIG. 3: D/S ratios in 61 lj-k and ai pTv decays, for ^Po 
and C^Po models. 



Acknowledgments 



technique, known as the Fock-Tani formahsm, long used 
in atomic physics. This formahsm has been apphed to 
hadron-hadron scattering interactions with constituent 
interchange. The chaUenge, resided in extending the ap- 
proach to include meson decay. After demonstrating that 
in lower order the result obtained was equivalent to the 
^Pq model, an additional feature pointed out was the ap- 
pearance of bound-state corrections in the effective decay 
Hamiltonian. These corrections present a natural modifi- 
cation in the qq interaction strength. As an example, we 
studied two decay processes 61 ojtt and ai pn. The 
D/S ratios, in Fig. ([3]), show that a common range of f3 
values for mesons is obtained. In a new calculation with 
the inclusion of other decay processes it might require 
different P values [2l|. The corrected model presents an 
interesting feature that for these two mesons the decay 
width differs slightly when compared with the '^Pq, but 
D/S ratios are improved. The examples studied here 
are encouraging but a more extensive survey of the light 
meson sector would be a necessary next step. The in- 
clusion of the full meson octet, in the evaluation of the 
bound-state correction, may provide a fine tuning for the 
model. 



The authors would like to thank H. Stocker, J. Aiche- 
lin and W. Greiner for important and enlightening discus- 
sions. This research was supported by Conselho Nacional 
de Desenvolvimento Cientifico e Tecnologico (CNPq), 
Universidade Federal do Rio Grande do Sul (UFRGS) 
and Universidade Federal de Pelotas (UFPel). 



APPENDIX A: SECOND AND THIRD ORDER 
OPERATORS 



The second order operators 
1 



where 



Map — rnl^Mp sin t cos t — mj^mp sin^ t 
-MlMp{2 - 2 cost - siii^ t) 
—M^mp (2 sin t — sin t cos t) . 



(AI) 



(A2) 
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The third order operators are 



where 



Act 



(B6) 



._$*pCT^^..^pCT.-t^-t^-^^^^ 



<^lP-<^>i^^-<^>P^^-qlqlqpMp, (A3) 



where 



Ma/Sj — ml^m/^my sin^ i + M^M^m^ (sini — sin'^ 
+MlmpMy (2 sin t - sin i cos t - sin^ i) 
+ [Ml^mpm-^ + rn^^Mpm^^ (— cost + cos'^ i) 
+mjjm/3M^ (— cos t + cos'^ i + sin^ t) 
+MlMpM^ (2 - cos i - cos^ t - sin^ t) 
+mJjM^M^ (sin i — sin t cos i — sin"^ t) 
= 2M;3(l-cosi) + m/3sint. (A4) 

APPENDIX B: THE SH^ HAMILTONIAN 

The SHm Hamihonian is evaluated from Eq. (|¥71) . The 
qj}'^^(0'^^ combination can be obtained from (jAip and 



= ql(^''>ql^^^V^,^Sfr{a,/3,j)V^,mimlm^ (Bl) 



with 



The qj}"^"^ qt''^'^ combination has an important feature: a 
contribution from a higher order operator. A new gen- 
erator Ma has to be evaluated, with the inclusion of the 
following fourth order term 

Mi^'> = lA^pAp^M^ - Im}^ [A„^, Aps] MsM^ 



M^Ms- 



The only relevant term in the q]}"^^ for meson decay is 

(B4) 

The resulting contribution is then 

5H2 - gt(i)^t(4) y^^^ ^ Sf^-'{a,P,-f)V^, mimlm^ (B5) 



(B3) 



The qjj^'^qt^^^ combination implies in a fifth order gen- 
erator to obtain the complete set of equations of motion 
(HZl) and (HH) 

Mt(5) = -Mlzl^M, + hllMlwl^^,MsM,iB7) 



where 



Z 



Q7/3 

a.Lu(3''f6 



Ia,,tL.-^-tL,a, 



[Ml 



A t'^ 



}ap-f5 — -2 [^07' A/3'5] 



57/3 



(B8) 



The only relevant terms in the g/^ for meson decay are 

-^A{pr-Xu;)^r<^r^''' 
X Aft(o) (t)Af^(°) (t)9.M^°) (t). (B9) 
The resulting contribution is 



(BIO) 



where 



The complete ^i?m Hamiltonian is 

SH^ = SH1 + SH2 + SH3 



Sr''{a,p,^)V,,,mimlm^ (B12) 



with 



Act 

7 



-<i>r*rA(pr;Ai.)$^^ 



-$r<i>;'"'A(pr;^A)$f . 



(B13) 
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APPENDIX C: DECAY, WAVE FUNCTION AND 
MATRIX ELEMENTS 

Wc will use the decay h'^{+z) uj{+z)Tr^ to illustrate 
the nature of our formalism and, simply quote the other 
case in the text. 



1. The wave function 

The general meson wave function can be written as 

a direct product of the spin Xs [the indexes s\ and 

,S2 arc the quark (antiquark) spin projections with (s = 
1 =>t and s = 2 ^i); the index Sa denotes the me- 
son spin]; flavor /J^^; color C^^'^ and space ^^f~^^~^^ 
components. 

In all our calculations the color component will be 
given by 



(JC1C2 _ 

\/3 



(C2) 



The spatial part is defined as harmonic oscillator wave 
functions 



where ^niiPi^Pj) is given by 
^nl{Pi,Pj) = Nni\pi- Pj\'' exp 



{Pi-Pjf 
8/32 



X Cn 



4/32 



i^(m(%,-p-) (C4) 



with the internal momentum, the spherical harmonic 

Yim, (3 a scale parameter, N^i the normalization constant 
dependent on the radial and orbital quantum numbers 



Nnl = 



2(n!) 



/33 r(n - 



3/2) 



(C5) 



The Laguerre polynomials ^ (p) are defined as 



^ (_)fer(n + Z + 3/2)("-'=) , 



c':hp)=y. 



fc=0 



k\T{k + l + Z/2) 



/ . (C6) 



In this paper two kinds of light non-strange mesons will 
be studied: 



1. Lqq = 



ip{p) = $00 (p) 



7^3/4^3/2 



exp 



'8/32 



(C7) 



2. Lgg = 1 



with 



2 


2 


p2 




p exp 





(C8) 

(C9) 



Returning to our example the pion, has J = and bi 
J =1. We choose the {+z) direction for this calculation, 
so the spin wave functions become 

= ^ (It I) -lit)) 

1'^) = |Tt) 

k) = ^ (It I) -lit)) (CIO) 

or in the x notation 

11 _ 1 12 _ 21 _ 22 _ n 

Xcj ' Xcj Xc^ Xcj ^ 

xilfel = x^!^. = ; xiX = -xl'M = ^- (cii) 

The flavor component /^^ " follows the same logic as the 
spin part 



|6+) = \7T+) = -\ud) 

\^) = -^{\u^) + \dd)) 



(C12) 



The bf and the 7r+ mesons have the same flavor contri- 
bution 

For oj, one has 

1 

71' 



J OJ J i 



2. The spin matrix elements 

In the evaluation of a decay amplitude, the following 
spin matrix element is necessary 



Xs 



(C15) 



with 



xi = I ) ; >^2 = I ^ I ; xi = I 1 I ; X5 - 1 



1 

I 

(CI6) 



By direct calculation one can show 

xl{s-p)xl = -Pz 

X*2{^-P)X2 = -{P.+iPy). (C17) 
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3. Matrix elements: 6^ lotv^ Decay 

The transition considered is of the form — > + 
nip, where the initial state is \A) = ml^ |0) and the final 

state is given by \BC) = ml^rnp |0). The matrix element 
of the uncorrected part results in 



{BC\H„,\A) ^-di-d2, 
di and ^2 are defined as 

di EE i;,, 
d2 = V. 



(C18) 



(C19) 



Equations (|C19p can be decomposed according to the 
sector of the wave function they correspond: flavor, color, 
spin-space: 



di = didldl- 
d2 = d^d^^d^-^ 



(C20) 



The matrix elements of the bound-state correction re- 
fer to diagrams (b), (c) and (d) of figure ([1]). The bound- 
state kernel's definition as <i>^'^$*'*"' implies in an addi- 
tional element, due to the contraction in the a index, a 
sum over species requirement jlJl . A question that nat- 
urally arises is, which states to include in this sum? We 
shall adopt in our calculation a restrictive choice: include 
in the sum only the particles that are present in the final 
state. For the 6^ decay, A (pr; Ai^) will have two con- 
tributions: w and 7r+. Similarly, the decay shall be 
corrected by the final state mesons p+ and 7r°. 

Due to the parity assignment of the spatial part, the in- 
tegration of diagram ^p) is zero. Spatial symmetry also 
implies that the matrix elements of diagrams ([TJ;) and 
(UJi) are equal. This simplifies our calculation, reducing 
the problem to the evaluation of diagram ([1]:) only. The 
bound-state correction (bsc) matrix element reduces to 
evaluate the following expression 



{BC\5Hra\A) 



if" 



jbsc 



(C21) 



4. — > uj-k'^ Decay (uncorrected) 
• Flavor: 



"l — "2 — Jit Ju Jbl "f^f.^ 

• Color: 



1 

V2- 



(C24) 



(C25) 



• Spin-space: 

The spin matrix element is 

(C26) 

where 

{Pl^vD - -lS{pl+p'u)xl [a • (p; -p;)]xf . (C27) 

Using (|C17p to evaluate (jC27p and after integrating mo- 
mentum conservation deltas one arrives in 

dl'^ = -7 / d^K (K^ - iKy) Lp(P-2K 



{C28) 

Introducing the spatial wave function and integrating 



-97/2 

"1 - OF./7. 



7 



35/2 y \^ 7^1/4/31/2 



1 - -x' 
9 



(C29) 



d| is obtained from d^ hy P ^ —P. The decay 
amplitude results 



7 



+ 



/ 2^ 
2 



1 2 2 
1 X 

9 



32^5 



X Y20 (Qx) f ei{x). 



(C30) 



where 



5. bf lutt'^ Decay (bound-state corrected) 



with 



dT - -^(dl^+dlrr) 
df" = J(d2.+rf2.) 



In (|C23p j refers to mesons u and tt^ 



(C22) 



d2j = K''^^J^A,{pT;X,y)<S>^''V^, = dy%d 



'S — e 
2j ■ 



(C23) 



The quantities between [. . .] in the following expres- 
sions are related to the bound-state kernel. 
• Flavor: 



/ 


^ fp fa Jt 

J TT J UJ 


f JpJt ffxffi' 
J UJ J UJ 


Jbl 


f 






ffxfa 

Jbl 


/ 

2u 


— ffp-f-r f fpftr 
J TT J UJ 


f SpSt tSxi^' 

J UJ J UJ 


ffxfcr 

Jbl 


/ 

27r 


— ffp-fT f fpfty 
J TT J UJ 




ffxfa 

Jbl 



2V2 
1 

1 
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Color: 



After integration one finds 



_ ^ SjCpC^ ^tc^Ct rirCpCT jtcac^i rcAC„ rc,,c„ 

1 



3V3' 
• Spin-space: 

The spin matrix element is 



(C32) 



dL = 
dL = 



Att Acj LAcj An; J A^i 



Att Aoj LAtt Att J Af);^ 



Att Ao; LAo; An; J Afaj^ s 



Att Acj LAtt Att J Abi s 



s — e 



s — e 



(C33) 



Due to symmetries in the spatial part the following rela- 
tions are true 



1 



d^Z' = ^. (C34) 



where dL is given by 

d\-' - I y" d^Kd^q{K,-iKy) 



. (2g) Fii (r!2,-) . 

(C35) 



7 



211/2 

75/2 y \^ 7^1/4 /31/2 

2^/2 n 



1- — 

21 



21 V 10 



x^Y2o i^x) } 62(2;). 



(C36) 



The decay amplitude for the bound-state correction 



;^bsc _ 



7 



27/2 > 

+ ^x^Y2o {n^)> e2{x). 



The total amplitude will be 



1 

21 



(C37) 



lC3P0 _ l j_o h -- — 



bsc _ 7 ^C3P0 



which are expressions ([50)) and (|5T|) . 



(C38) 



[1] F. E. Close and A. Kirk, Eur. Phys. J. C 21, 531 (2001). 
[2] U. Thoma, Eur. Phys. J. A 18, 135 (2003). [11] 
[3] E. Klempt, A. Zaitsev, Phys. Rept. 454, 1 (2007). 
[4] T. Barnes, Int. J. Mod. Phys. A 21, 5583 (2006). 
[5] L. Micu, Nucl. Phys. B 10, 512 (1969). [12] 
[6] A. Leyaouanc, L. Oliver, O. Pene and J. Raynal, Phys. 
Rev. D 8, 2223 (1973); tbtd. D 9, 1415 (1974); D 11, 680 
(1975); D 11, 1272 (1975); Phys. Rev. Lett. B 71, 397 [13] 
(1977). 

[7] E. S. Ackleh, T. Barnes and E. S. Swanson, Phys. Rev. 

D 54, 6811 (1996). [14] 
[8] T. Barnes, F. E. Close, P. R. Page and E. S. Swanson, 

Phys. Rev. D 55, 4157 (1997). [15] 
[9] T. Barnes, N. Black and P.R. Page, Phys. Rev. D 68, [16] 

054014 (2003). 

[10] T. Barnes, S. Godfrey and E. S. Swanson, Phys. Rev. D [17] 



72, 054026 (2005). 

M.D. Girardeau, Phys. Rev. Lett. 27, 1416 (1971), ibid. 
J. Math. Phys. 16, 1901 (1975); Phys Rev. A 26, 217 
(1982). 

D. Hadjimichef, G. Krein, S. Szpigel and J. S. da Veiga, 
Ann. of Phys. 268, 105 (1998); ibid. Phys. Lett. B 367, 
317 (1996). 

S. Szpigel, Interagao Meson-Meson no Formalismo Fock- 
Tani. PhD thesis, Instituto de Fisica, Universidade de 
Sao Paulo, Sao Paulo, 1995. 

D. Hadjimichef, J. Haidenbauer, G. Krein, Phys. Rev. C 
63, 035204 (2001); ibid. C 66, 055214 (2002). 

D. T. da Silva, D. Hadjimichef, J. Phys. G 30, 191 (2004). 
M. L. L. Silva, D. Hadjimichef, C. A. Z. Vasconcellos, B. 

E. J. Bodmann, J. Phys. G 32, 475 (2006). 

S. Weinberg, Phys. Rev. 130, 776 (1963); ibid. 131, 440 



12 



(1963); M. Scadron and S. Weinberg, Phys. Rev. 133, 

B1589 (1964); M. Scadron, S. Weinberg and J. Wright, 

Phys. Rev. 135 B202 (1964). 
[18] M. Oka and K. Yazaki, Prog. Theor. Phys. 66 556 

(1981); ibid. 572 (1981). 
[19] E. S. Swanson, Ann. of Phys. 220, 73 (1992); T. Barnes 

and E. S. Swanson, Phys. Rev. D 46, 131 (1992); T. 



Barnes, S. Capstick, M. D. Kovarik and E. S. Swanson, 
Phys. Rev. C 48, 539 (1993). 
[20] Particle Data Group, W.-M. Yao et al, J. Phys. G 33, 1 
(2006). 

[21] M. Strohmeier-Presicek, T. Gutsche, R. Vinh Mau, and 
A. Faessler, Phys. Rev. D 60, 054010 (1999). 



